Transcriptomic Analysis Reveals Sixteen Potential Genes Associated with the Successful Differentiation of Antibody-Secreting Cells through the Utilization of Unfolded Protein Response Mechanisms in Robust Responders to the Influenza Vaccine.

The seasonal influenza vaccine remains one of the vital recommended infection control measures for the elderly with chronic illnesses. We investigated the immunogenicity of a single dose of influenza vaccine in 123 seronegative participants and classified them into four distinct groups, determined by the promptness of vaccine response, the longevity of humoral immunity, and the likelihood of exhibiting cross-reactivity. Subsequently, we used transcriptional profiling and differential gene expression analysis to identify potential genes directly associated with the robust response to the vaccine. The group of exemplary vaccine responders differentially expressed 16 genes, namely: MZB1, MYDGF, TXNDC5, TXNDC11, HSP90B1, FKBP11, PDIA5, PRDX4, CD38, SDC1, TNFRSF17, TNFRSF13B, PAX5, POU2AF1, IRF4, and XBP1. Our findings point out a list of expressed proteins that are related to B cell proliferation, unfolded protein response, and cellular haemostasis, as well as a linkage of these expressions to the survival of long-lived plasma cells.

Genetic architecture of alcohol consumption identified by a genotype-stratified GWAS and impact on esophageal cancer risk in Japanese people.

An East Asian-specific variant on aldehyde dehydrogenase 2 (ALDH2 rs671, G>A) is the major genetic determinant of alcohol consumption. We performed an rs671 genotype-stratified genome-wide association study meta-analysis of alcohol consumption in 175,672 Japanese individuals to explore gene-gene interactions with rs671 behind drinking behavior. The analysis identified three genome-wide significant loci (GCKR, KLB, and ADH1B) in wild-type homozygotes and six (GCKR, ADH1B, ALDH1B1, ALDH1A1, ALDH2, and GOT2) in heterozygotes, with five showing genome-wide significant interaction with rs671. Genetic correlation analyses revealed ancestry-specific genetic architecture in heterozygotes. Of the discovered loci, four (GCKR, ADH1B, ALDH1A1, and ALDH2) were suggested to interact with rs671 in the risk of esophageal cancer, a representative alcohol-related disease. Our results identify the genotype-specific genetic architecture of alcohol consumption and reveal its potential impact on alcohol-related disease risk.

Trivalent inactivated influenza vaccine response and immunogenicity assessment after one week and three months in repeatedly vaccinated adults.

BACKGROUND: The influenza vaccine administrated every year is a recommended infection control procedure for individuals above the age of six months. However, the effectiveness of repeated annual vaccination is still an active research topic. Therefore, we investigated the vaccine immunogenicity in two independent groups: previously vaccinated versus non-vaccinated individuals at three time points; prior vaccination, one week and three months post vaccination. The assessment enabled us to evaluate the elicited immune responses and the durability of the induced protection in both groups. RESEARCH DESIGN AND METHODS: A research study was conducted to assess the immunogenicity of a single dose of Trivalent Inactivated Influenza Vaccine (A/H1N1, A/H3N2, and B) in 278 healthy adults aged between 32 and 66 years. Almost half of the participants, 140 (50·36%), received influenza vaccination at least once precursor to past influenza seasons. One blood sample was taken prior to vaccination for complete blood analysis and baseline immunogenicity assessment. The selected study participants received a single vaccine dose on the first day, and then followed up for three months. Two blood samples were taken after one week and three months post vaccination, respectively, for vaccine immunogenicity assessment. RESULTS: Before vaccination, the seroprotection, defined as a hemagglutination-inhibiting titer of =>1:40, was detected for the three vaccine virus strains in 20 previously vaccinated participants (14·29%) [8·95%, 21·2%]. We compared the overall vaccine response for the three virus strains using a normalized response score calculated from linearly transformed titer measurements; the score before vaccination was 84% higher in the previously vaccinated group and the mean difference between the two groups was statistically significant. Three months post-vaccination, we didn't find a significant difference in vaccine responses; the number of fully seroprotected individuals became 48 (34·29%) [26·48%, 42·77%] in the previously vaccinated group and 59 (42·75%) [34·37%, 51·45%] in the non-vaccinated group. The calculated response score was almost equal in both groups and the mean difference was no longer statistically significant. CONCLUSION: Our findings suggest that a single dose of influenza vaccine is equally protective after three months for annually vaccinated adults and first-time vaccine receivers.

Improved genetic prediction of the risk of knee osteoarthritis using the risk factor-based polygenic score.

BACKGROUND: Polygenic risk score (PRS) analysis is used to predict disease risk. Although PRS has been shown to have great potential in improving clinical care, PRS accuracy assessment has been mainly focused on European ancestry. This study aimed to develop an accurate genetic risk score for knee osteoarthritis (OA) using a multi-population PRS and leveraging a multi-trait PRS in the Japanese population. METHODS: We calculated PRS using PRS-CS-auto, derived from genome-wide association study (GWAS) summary statistics for knee OA in the Japanese population (same ancestry) and multi-population. We further identified risk factor traits for which PRS could predict knee OA and subsequently developed an integrated PRS based on multi-trait analysis of GWAS (MTAG), including genetically correlated risk traits. PRS performance was evaluated in participants of the Nagahama cohort study who underwent radiographic evaluation of the knees (n = 3,279). PRSs were incorporated into knee OA integrated risk models along with clinical risk factors. RESULTS: A total of 2,852 genotyped individuals were included in the PRS analysis. The PRS based on Japanese knee OA GWAS was not associated with knee OA (p = 0.228). In contrast, PRS based on multi-population knee OA GWAS showed a significant association with knee OA (p = 6.7 × 10-5, odds ratio (OR) per standard deviation = 1.19), whereas PRS based on MTAG of multi-population knee OA, along with risk factor traits such as body mass index GWAS, displayed an even stronger association with knee OA (p = 5.4 × 10-7, OR = 1.24). Incorporating this PRS into traditional risk factors improved the predictive ability of knee OA (area under the curve, 74.4% to 74.7%; p = 0.029). CONCLUSIONS: This study showed that multi-trait PRS based on MTAG, combined with traditional risk factors, and using large sample size multi-population GWAS, significantly improved predictive accuracy for knee OA in the Japanese population, even when the sample size of GWAS of the same ancestry was small. To the best of our knowledge, this is the first study to show a statistically significant association between the PRS and knee OA in a non-European population. TRIAL REGISTRATION: No. C278.

Identification of risk factors for elevated serum IgG4 levels in subjects in a large-scale health checkup cohort study.

Introduction: To allow the identification of IgG4-related disease (IgG4-RD) from a subclinical phase as it is important to understand the risk of elevated serum IgG4 levels. We planned to evaluate serum IgG4 levels in the participants of the Nagasaki Islands Study (NaIS), a large-scale health checkup cohort study. Methods: This study included 3,240 individuals who participated in the NaIS between 2016 and 2018 and consented to participate in the study. Serum IgG4, IgG, and IgE levels and human leukocyte antigen (HLA) genotyping results of the NaIS subjects as well as lifestyle habits and peripheral blood test results were analyzed. The magnetic bead panel assay (MBA) and the standard nephelometry immunoassay (NIA) were used to measure serum IgG4 levels. The data were evaluated using multivariate analysis to identify lifestyle and genetic factors associated with elevated serum IgG4 levels. Results: Serum IgG4 levels measured with the NIA and MBA showed a tight positive correlation between the two groups (correlation coefficient 0.942). The median age of the participants in the NaIS was 69 years [63-77]. The median serum IgG4 level was 30.2 mg/dL [IQR 12.5-59.8]. Overall, 1019 (32.1%) patients had a history of smoking. When the subjects were stratified into three groups based on the smoking intensity (pack-year), the serum IgG4 level was significantly higher among those with a higher smoking intensity. Accordingly, the multivariate analysis identified a significant relationship between smoking status and serum IgG4 elevation. Conclusion: In this study, smoking was identified as a lifestyle factor correlating positively with elevated serum IgG4 levels.

Genetic influences on human blood metabolites in the Japanese population.

An increase in ethnic diversity in genetic studies has the potential to provide unprecedented insights into how genetic variations influence human phenotypes. In this study, we conducted a quantitative trait locus (QTL) analysis of 121 metabolites measured using gas chromatography-mass spectrometry with plasma samples from 4,888 Japanese individuals. We found 60 metabolite-gene associations, of which 13 have not been previously reported. Meta-analyses with another Japanese and a European study identified six and two additional unreported loci, respectively. Genetic variants influencing metabolite levels were more enriched in protein-coding regions than in the regulatory regions while being associated with the risk of various diseases. Finally, we identified a signature of strong negative selection for uric acid ( S ˆ  = -1.53, p = 6.2 × 10-18). Our study expanded the knowledge of genetic influences on human blood metabolites, providing valuable insights into their physiological, pathological, and selective properties.

Genome-Wide Association Study of Intracranial Artery Stenosis Followed by Phenome-Wide Association Study.

The genetic background of intracranial artery stenosis (ICAS), a major cause of ischemic stroke, remains elusive. We performed the world's first genome-wide association study (GWAS) of ICAS using DNA samples from Japanese subjects, to identify the genetic factors associated with ICAS and their correlation with clinical features. We also conducted a phenome-wide association study (PheWAS) of the top variant identified via GWAS to determine its association with systemic disease. The GWAS involved 408 patients with ICAS and 349 healthy controls and utilized an Asian Screening Array of venous blood samples. The PheWAS was performed using genotypic and phenotypic data of the Biobank Japan Project, which contained information on 46 diseases and 60 quantitative trait data from > 150,000 Japanese individuals. The GWAS revealed that the East Asian-specific functional variant of RNF213, rs112735431 (c.14429G > A, p.Arg4810Lys), was associated with ICAS (odds ratio, 12.3; 95% CI 5.5 to 27.5; P = 7.8 × 10-10). Stratified analysis within ICAS cases demonstrated that clinical features of those with and without the risk allele were different. PheWAS indicated that high blood pressure and angina were significantly associated with RNF213 rs112735431. The first GWAS of ICAS, which stratifies subpopulations within the ICAS cases with distinct clinical features, revealed that RNF213 rs112735431 was the most significant variant associated with ICAS. Thus, RNF213 rs112735431 shows potential as an important clinical biomarker that characterizes pleiotropic risk in various vascular diseases, such as blood pressure and angina, thereby facilitating personalized medicine for systemic vascular diseases in East Asian populations.

Multi-ancestry genome-wide association analyses identify novel genetic mechanisms in rheumatoid arthritis.

Rheumatoid arthritis (RA) is a highly heritable complex disease with unknown etiology. Multi-ancestry genetic research of RA promises to improve power to detect genetic signals, fine-mapping resolution and performances of polygenic risk scores (PRS). Here, we present a large-scale genome-wide association study (GWAS) of RA, which includes 276,020 samples from five ancestral groups. We conducted a multi-ancestry meta-analysis and identified 124 loci (P < 5 × 10-8), of which 34 are novel. Candidate genes at the novel loci suggest essential roles of the immune system (for example, TNIP2 and TNFRSF11A) and joint tissues (for example, WISP1) in RA etiology. Multi-ancestry fine-mapping identified putatively causal variants with biological insights (for example, LEF1). Moreover, PRS based on multi-ancestry GWAS outperformed PRS based on single-ancestry GWAS and had comparable performance between populations of European and East Asian ancestries. Our study provides several insights into the etiology of RA and improves the genetic predictability of RA.

A common deletion at BAK1 reduces enhancer activity and confers risk of intracranial germ cell tumors.

Intracranial germ cell tumors (IGCTs) are rare brain neoplasms that mainly occur in children and adolescents with a particularly high incidence in East Asian populations. Here, we conduct a genome-wide association study (GWAS) of 133 patients with IGCTs and 762 controls of Japanese ancestry. A common 4-bp deletion polymorphism in an enhancer adjacent to BAK1 is significantly associated with the disease risk (rs3831846; P = 2.4 × 10-9, odds ratio = 2.46 [95% CI: 1.83-3.31], minor allele frequency = 0.43). Rs3831846 is in strong linkage disequilibrium with a testicular GCTs susceptibility variant rs210138. In-vitro reporter assays reveal rs3831846 to be a functional variant attenuating the enhancer activity, suggesting its contribution to IGCTs predisposition through altering BAK1 expression. Risk alleles of testicular GCTs derived from the European GWAS show significant positive correlations in the effect sizes with the Japanese IGCTs GWAS (P = 1.3 × 10-4, Spearman's ρ = 0.48). These results suggest the shared genetic susceptibility of GCTs beyond ethnicity and primary sites.

The association between the Moyamoya disease susceptible gene RNF213 variant and incident cardiovascular disease in a general population: the Nagahama study.

OBJECTIVE: An association between the Moyamoya disease susceptible gene ring finger protein 213 (RNF213) variant and ischemic stroke and coronary artery disease has been suggested in case-control studies. We aimed to investigate the possible association between the RNF213 variant and the incidence of cardiovascular disease in a general population. METHODS: The study participants consisted of 9153 Japanese community residents without history of cardiovascular disease. The clinical parameters employed in this analysis were observed at baseline between 2008 and 2010. The RNF213 p.R4859K variant was determined by TaqMan probe assay and then confirmed by Sanger sequencing. RESULTS: During 8.52 years follow-up period, we observed 214 incident cases of cardiovascular diseases (99 total stroke cases, 119 major adverse cardiac event cases, including 4 cases of both). The incidence rate was higher for the variant allele carriers (120 cases; incidence rate, 71.0 per 10 000 person-years) than for the homozygotes of the wild-type allele (26.9), and the group differences achieved statistical significance (P = 0.009). Although the RNF213 variant was also associated with systolic blood pressure (dominant model: coefficient of 8.19 mmHg; P < 0.001), the Cox regression analysis adjusted for major covariates including systolic blood pressure identified the RNF213 variant as an independent determinant for cardiovascular disease (hazard ratio of 3.41, P = 0.002) and major adverse cardiac event (hazard ratio of 3.80, P = 0.010) but not with total stroke (P = 0.102). CONCLUSION: The Moyamoya disease susceptible RNF213 variant was associated with blood pressure and the incidence of cardiovascular disease in a Japanese general population.

Genome wide association study of HTLV-1-associated myelopathy/tropical spastic paraparesis in the Japanese population.

HTLV-1-associated myelopathy (HAM/TSP) is a chronic and progressive inflammatory disease of the central nervous system. The aim of our study was to identify genetic determinants related to the onset of HAM/TSP in the Japanese population. We conducted a genome-wide association study comprising 753 HAM/TSP patients and 899 asymptomatic HTLV-1 carriers. We also performed comprehensive genotyping of HLA-A, -B, -C, -DPB1, -DQB1, and -DRB1 genes using next-generation sequencing technology for 651 HAM/TSP patients and 804 carriers. A strong association was observed in HLA class I (P = 1.54 × 10-9) and class II (P = 1.21 × 10-8) loci with HAM/TSP. Association analysis using HLA genotyping results showed that HLA-C*07:02 (P = 2.61 × 10-5), HLA-B*07:02 (P = 4.97 × 10-10), HLA-DRB1*01:01 (P = 1.15 × 10-9) and HLA-DQB1*05:01 (P = 2.30 × 10-9) were associated with disease risk, while HLA-B*40:06 (P = 3.03 × 10-5), HLA-DRB1*15:01 (P = 1.06 × 10-5) and HLA-DQB1*06:02 (P = 1.78 × 10-6) worked protectively. Logistic regression analysis identified amino acid position 7 in the G-BETA domain of HLA-DRB1 as strongly associated with HAM/TSP (P = 9.52 × 10-10); individuals homozygous for leucine had an associated increased risk of HAM/TSP (odds ratio, 9.57), and proline was protective (odds ratio, 0.65). Both associations were independent of the known risk associated with proviral load. DRB1-GB-7-Leu was not significantly associated with proviral load. We have identified DRB1-GB-7-Leu as a genetic risk factor for HAM/TSP development independent of proviral load. This suggests that the amino acid residue may serve as a specific marker to identify the risk of HAM/TSP even without knowledge of proviral load. In light of its allele frequency worldwide, this biomarker will likely prove useful in HTLV-1 endemic areas across the globe.

Genetic determinants of risk in autoimmune pulmonary alveolar proteinosis.

Pulmonary alveolar proteinosis (PAP) is a devastating lung disease caused by abnormal surfactant homeostasis, with a prevalence of 6-7 cases per million population worldwide. While mutations causing hereditary PAP have been reported, the genetic basis contributing to autoimmune PAP (aPAP) has not been thoroughly investigated. Here, we conducted a genome-wide association study of aPAP in 198 patients and 395 control participants of Japanese ancestry. The common genetic variant, rs138024423 at 6p21, in the major-histocompatibility-complex (MHC) region was significantly associated with disease risk (Odds ratio [OR] = 5.2; P = 2.4 × 10-12). HLA fine-mapping revealed that the common HLA class II allele, HLA-DRB1*08:03, strongly drove this signal (OR = 4.8; P = 4.8 × 10-12), followed by an additional independent risk allele at HLA-DPß1 amino acid position 8 (OR = 0.28; P = 3.4 × 10-7). HLA-DRB1*08:03 was also associated with an increased level of anti-GM-CSF antibody, a key driver of the disease (ß = 0.32; P = 0.035). Our study demonstrated a heritable component of aPAP, suggesting an underlying genetic predisposition toward an abnormal antibody production.

Genome-wide association meta-analysis identifies GP2 gene risk variants for pancreatic cancer.

Pancreatic cancer is the fourth leading cause of cancer-related deaths in Japan. To identify risk loci, we perform a meta-analysis of three genome-wide association studies comprising 2,039 pancreatic cancer patients and 32,592 controls in the Japanese population. Here, we identify 3 (13q12.2, 13q22.1, and 16p12.3) genome-wide significant loci (P < 5.0 × 10-8), of which 16p12.3 has not been reported in the Western population. The lead single nucleotide polymorphism (SNP) at 16p12.3 is rs78193826 (odds ratio = 1.46, 95% confidence interval = 1.29-1.66, P = 4.28 × 10-9), an Asian-specific, nonsynonymous glycoprotein 2 (GP2) gene variant. Associations between selected GP2 gene variants and pancreatic cancer are replicated in 10,822 additional cases and controls of East Asian origin. Functional analyses using cell lines provide supporting evidence of the effect of rs78193826 on KRAS activity. These findings suggest that GP2 gene variants are probably associated with pancreatic cancer susceptibility in populations of East Asian ancestry.

Association of ALPL variants with serum alkaline phosphatase and bone traits in the general Japanese population: The Nagahama Study.

Although alkaline phosphatase (ALP) activity is relatively low in carriers of recessive type hypophosphatasia (HPP), most are asymptomatic and therefore do not undergo medical evaluations. We analyzed the association of ALP-encoding ALPL variants with serum ALP and bone traits in the general Japanese population. Study participants (n = 9671) were from the Nagahama Study, which was a longitudinal cohort study of an apparently healthy general Japanese population. ALPL variants were analyzed by whole-genome sequencing or TaqMan probe assays using DNA extracted from peripheral blood samples. The speed of sound in calcaneal bone was assessed by quantitative ultrasound (QUS) and used as surrogate measures of bone mineral density. We identified 13 ALPL variants. Minor allele frequencies of three variants were higher than expected. Variant c.529G > A has been reported as a possible pathogenic variant for adult type HPP. Variants c.979C > T and c.1559delT are reported as pathogenic variants for perinatal severe HPP or infantile HPP. The allele frequencies of c.529G > A, c.979C > T, and c.1559delT were 0.0107, 0.0040, and 0.0014, respectively. Serum ALP activity was significantly lower and differed among the three variants (P < 0.001), as well as between individuals with and without any of the three variants (P < 0.001). Serum ALP activity was inversely associated with QUS values, although no direct association was observed between the ALPL variants and QUS values. An association between serum ALP activity and QUS was confirmed; however, we failed to detect an association between ALPL variants and bone traits in the general Japanese population.

Genetic basis for plasma amino acid concentrations based on absolute quantification: a genome-wide association study in the Japanese population.

To assess the use of plasma free amino acids (PFAAs) as biomarkers for metabolic disorders, it is essential to identify genetic factors that influence PFAA concentrations. PFAA concentrations were absolutely quantified by liquid chromatography-mass spectrometry using plasma samples from 1338 Japanese individuals, and genome-wide quantitative trait locus (QTL) analysis was performed for the concentrations of 21 PFAAs. We next conducted a conditional QTL analysis using the concentration of each PFAA adjusted by the other 20 PFAAs as covariates to elucidate genetic determinants that influence PFAA concentrations. We identified eight genes that showed a significant association with PFAA concentrations, of which two, SLC7A2 and PKD1L2, were identified. SLC7A2 was associated with the plasma levels of arginine and ornithine, and PKD1L2 with the level of glycine. The significant associations of these two genes were revealed in the conditional QTL analysis, but a significant association between serine and the CPS1 gene disappeared when glycine was used as a covariate. We demonstrated that conditional QTL analysis is useful for determining the metabolic pathways predominantly used for PFAA metabolism. Our findings will help elucidate the physiological roles of genetic components that control the metabolism of amino acids.

Genetic determinants and an epistasis of LILRA3 and HLA-B*52 in Takayasu arteritis.

Takayasu arteritis (TAK) is a systemic vasculitis with severe complications that affects the aorta and its large branches. HLA-B*52 is an established susceptibility locus to TAK. To date, there are still only a limited number of reports concerning non-HLA susceptibility loci to TAK. We conducted a genome-wide association study (GWAS) and a follow-up study in a total of 633 TAK cases and 5,928 controls. A total of 510,879 SNPs were genotyped, and 5,875,450 SNPs were imputed together with HLA-B*52. Functional annotation of significant loci, enhancer enrichment, and pathway analyses were conducted. We identified four unreported significant loci, namely rs2322599, rs103294, rs17133698, and rs1713450, in PTK2B, LILRA3/LILRB2, DUSP22, and KLHL33, respectively. Two additional significant loci unreported in non-European GWAS were identified, namely HSPA6/FCGR3A and chr21q.22. We found that a single variant associated with the expression of MICB, a ligand for natural killer (NK) cell receptor, could explain the entire association with the HLA-B region. Rs2322599 is strongly associated with the expression of PTK2B Rs103294 risk allele in LILRA3/LILRB2 is known to be a tagging SNP for the deletion of LILRA3, a soluble receptor of HLA class I molecules. We found a significant epistasis effect between HLA-B*52 and rs103294 (P = 1.2 × 10-3). Enhancer enrichment analysis and pathway analysis suggested the involvement of NK cells (P = 8.8 × 10-5, enhancer enrichment). In conclusion, four unreported TAK susceptibility loci and an epistasis effect between LILRA3 and HLA-B*52 were identified. HLA and non-HLA regions suggested a critical role for NK cells in TAK.

Prediction model for pancreatic cancer risk in the general Japanese population.

Genome-wide association studies (GWASs) have identified many single nucleotide polymorphisms (SNPs) that are significantly associated with pancreatic cancer susceptibility. We sought to replicate the associations of 61 GWAS-identified SNPs at 42 loci with pancreatic cancer in Japanese and to develop a risk model for the identification of individuals at high risk for pancreatic cancer development in the general Japanese population. The model was based on data including directly determined or imputed SNP genotypes for 664 pancreatic cancer case and 664 age- and sex-matched control subjects. Stepwise logistic regression uncovered five GWAS-identified SNPs at five loci that also showed significant associations in our case-control cohort. These five SNPs were included in the risk model and also applied to calculation of the polygenic risk score (PRS). The area under the curve determined with the leave-one-out cross-validation method was 0.63 (95% confidence interval, 0.60-0.66) or 0.61 (0.58-0.64) for versions of the model that did or did not include cigarette smoking and family history of pancreatic cancer in addition to the five SNPs, respectively. Individuals in the lowest and highest quintiles for the PRS had odds ratios of 0.62 (0.42-0.91) and 1.98 (1.42-2.76), respectively, for pancreatic cancer development compared with those in the middle quintile. We have thus developed a risk model for pancreatic cancer that showed moderately good discriminatory ability with regard to differentiation of pancreatic cancer patients from control individuals. Our findings suggest the potential utility of a risk model that incorporates replicated GWAS-identified SNPs and established demographic or environmental factors for the identification of individuals at increased risk for pancreatic cancer development.

Risk estimation model for nonalcoholic fatty liver disease in the Japanese using multiple genetic markers.

The genetic factors affecting the natural history of nonalcoholic fatty liver disease (NAFLD), including the development of nonalcoholic steatohepatitis (NASH) and NASH-derived hepatocellular carcinoma (NASH-HCC), are still unknown. In the current study, we sought to identify genetic factors related to the development of NAFLD, NASH, and NASH-HCC, and to establish risk-estimation models for them. For these purposes, 936 histologically proven NAFLD patients were recruited, and genome-wide association (GWA) studies were conducted for 902, including 476 NASH and 58 NASH-HCC patients, against 7,672 general-population controls. Risk estimations for NAFLD and NASH were then performed using the SNPs identified as having significant associations in the GWA studies. We found that rs2896019 in PNPLA3 [p = 2.3x10-31, OR (95%CI) = 1.85 (1.67-2.05)], rs1260326 in GCKR [p = 9.6x10-10, OR (95%CI) = 1.38(1.25-1.53)], and rs4808199 in GATAD2A [p = 2.3x10-8, OR (95%CI) = 1.37 (1.23-1.53)] were significantly associated with NAFLD. Notably, the number of risk alleles in PNPLA3 and GATAD2A was much higher in Matteoni type 4 (NASH) patients than in type 1, type 2, and type 3 NAFLD patients. In addition, we newly identified rs17007417 in DYSF [p = 5.2x10-7, OR (95%CI) = 2.74 (1.84-4.06)] as a SNP associated with NASH-HCC. Rs641738 in TMC4, which showed association with NAFLD in patients of European descent, was not replicated in our study (p = 0.73), although the complicated LD pattern in the region suggests the necessity for further investigation. The genetic variants of PNPLA3, GCKR, and GATAD2A were then used to estimate the risk for NAFLD. The obtained Polygenic Risk Scores showed that the risk for NAFLD increased with the accumulation of risk alleles [AUC (95%CI) = 0.65 (0.63-0.67)]. CONCLUSIONS: We demonstrated that NASH is genetically and clinically different from the other NAFLD subgroups. We also established risk-estimation models for NAFLD and NASH using multiple genetic markers. These models can be used to improve the accuracy of NAFLD diagnosis and to guide treatment decisions for patients.

Transethnic meta-analysis identifies GSDMA and PRDM1 as susceptibility genes to systemic sclerosis.

OBJECTIVES: Systemic sclerosis (SSc) is an autoimmune disease characterised by skin and systemic fibrosis culminating in organ damage. Previous genetic studies including genome-wide association studies (GWAS) have identified 12 susceptibility loci satisfying genome-wide significance. Transethnic meta-analyses have successfully expanded the list of susceptibility genes and deepened biological insights for other autoimmune diseases. METHODS: We performed transethnic meta-analysis of GWAS in the Japanese and European populations, followed by a two-staged replication study comprising a total of 4436 cases and 14 751 controls. Associations between significant single nuclear polymorphisms (SNPs) and neighbouring genes were evaluated. Enrichment analysis of H3K4Me3, a representative histone mark for active promoter was conducted with an expanded list of SSc susceptibility genes. RESULTS: We identified two significant SNP in two loci, GSDMA and PRDM1, both of which are related to immune functions and associated with other autoimmune diseases (p=1.4×10-10 and 6.6×10-10, respectively). GSDMA also showed a significant association with limited cutaneous SSc. We also replicated the associations of previously reported loci including a non-GWAS locus, TNFAIP3. PRDM1 encodes BLIMP1, a transcription factor regulating T-cell proliferation and plasma cell differentiation. The top SNP in GSDMA was a missense variant and correlated with gene expression of neighbouring genes, and this could explain the association in this locus. We found different human leukocyte antigen (HLA) association patterns between the two populations. Enrichment analysis suggested the importance of CD4-naïve primary T cell. CONCLUSIONS: GSDMA and PRDM1 are associated with SSc. These findings provide enhanced insight into the genetic and biological basis of SSc.

Genotype Analyses in the Japanese and Belarusian Populations Reveal Independent Effects of rs965513 and rs1867277 but Do Not Support the Role of FOXE1 Polyalanine Tract Length in Conferring Risk for Papillary Thyroid Carcinoma.

BACKGROUND: Several functional single-nucleotide polymorphisms (SNPs) at the FOXE1 locus on chromosome 9q22.33 have been associated with the risk for papillary thyroid carcinoma (PTC). This study set out to elucidate whether their effects are independent, using genotyping results in populations of Asian and European descent. METHODS: SNPs rs965513 and rs1867277 and a polymorphic region determining the length of the FOXE1 polyalanine (poly-Ala) tract were genotyped in 501 patients with PTC and 748 healthy individuals from Japan, and in 660 patients and 820 population controls from Belarus. Functional analysis of transactivation activities of FOXE1 isoforms with varying number of alanine repeats was performed by a Dual-Luciferase® Assay. RESULTS: All three polymorphisms were significantly associated with PTC in both populations on univariate analysis. However, conditional analysis revealed independent effects of rs965513 and rs1867277 SNPs but not of the FOXE1 poly-Ala polymorphism. The independent effect of the lead rs965513 SNP was observed in both populations, while that of rs1867277 was only identified in the Japanese population, in which linkage disequilibrium between the three polymorphisms is markedly weaker. Despite the strong decrease in transcriptional activity with increasing FOXE1 poly-Ala tract length, no difference in transactivation potential of the FOXE1 poly-Ala isoforms could be seen after adjustment for the minimal promoter activity in the reporter vectors. Plasmids encoding FOXE1 isoforms of increasing poly-Ala tract length were also found to produce less FOXE1 protein after cell transfection. CONCLUSIONS: The functional variants rs965513 and rs1867277 independently contribute to genetic predisposition to PTC, while a contributing role of the FOXE1 poly-Ala polymorphism could not be confirmed.